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Abstract 

We present a comprehensive analysis of coherence properties of the radiation from X- 
ray free electron laser (XFEL). We consider practically important case when XFEL is op- 
timized for maximum gain. Such an optimization allows to reduce significantly parameter 
space. Application of similarity techniques to the results of numerical simulations allows to 
present all output characteristics of the optimized XFEL as functions of the only parameter, 
ratio of the emittance to the radiation wavelength, e = 27re/A. Our studies show that op- 
timum performance of the XFEL in terms of transverse coherence is achieved at the value 
of the parameter e of about unity. At smaller values of e the degree of transverse coherence 
is reduced due to strong influence of poor longitudinal coherence on a transverse one. At 
large values of the emittance the degree of transverse coherence degrades due to poor mode 
selection. Comparative analysis of existing XFEL projects, European XFEL, LCLS, and 
SCSS is presented as well. 
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1 Introduction 



Free electron lasing at wavelengths shorter than ultraviolet can be achieved with a single- 
pass, high-gain FEL amplifier. Because a lack of powerful, coherent seeding sources short- 
wavelength FEL amplifiers work in so called Self-Amplified Spontaneous Emission (SASE) 
mode when amplification process starts from shot noise in the electron beam [1,2,3]. Present 
level of accelerator and FEL techniques holds potential for SASE FELs to generate wavelengths 
as short as 0.1 nm [4,5,6]. 

Experimental realization of X-ray FELs (XFELs) developed very rapidly during last decade. 
The first demonstration of the SASE FEL mechanism took place in 1997 in the infrared wave- 
length range [7]. In September of 2000, a group at Argonne National Laboratory (ANL) became 
the first to demonstrate saturation in a visible (390 nm) SASE FEL [8]. In September 2001, a 
group at DESY (Hamburg, Germany) has demonstrated lasing to saturation at 98 nm [9,10]. 
In June 2006 saturation has been achieved at 13 nm, the shortest wavelength ever generated 
by FELs. The experimental results have been achieved at FLASH ("F"ree-Electron-"LAS"er 
in "H"amburg). Regular user operation of FLASH started in 2005 [11]. Currently FLASH pro- 
duces GW- level, laser-like VUV radiation pulses with 10 to 50 fs duration in the wavelength 
range 13-45 nm. After the energy upgrade of the FLASH linac to 1 GeV planned in 2007, it will 
be possible to generate wavelengths down to 6 nm. 

Recently the German government, encouraged by these results, approved funding a hard X-ray 
SASE FEL user facility - the European X-Ray Free Electron Laser [4]. The US Department of 
Energy (DOE) has given SLAG the goahead for the engineering design of the Linac Coherent 
Light Source (LCLS) to be constructed at SLAG [5]. These devices should produce 100 fs X-ray 
pulses with over 10 to 100 GW of peak power. The main difference between projects is the linear 
accelerator, an existing room temperature linac for LGLS at SLAG, and future superconducting 
linac for the European XFEL. The XFEL based on superconducting accelerator technology will 
make possible not only a jump in a peak brilliance by ten orders of magnitude, but also increase 
by five orders of magnitude in average brilliance. The LGLS and European XFEL projects are 
scheduled to start operation in 2009 and 2013, respectively. 

In the X-ray FEL the radiation is produced by the electron beam during single-pass of the un- 
dulator [1,2,3]. The amplification process starts from the shot noise in the electron beam. Any 
random fluctuations in the beam current correspond to a modulation of the beam current at all 
frequencies simultaneously. When the electron beam enters the undulator, the presence of the 
beam modulation at frequencies close to the resonance frequency initiates the process of radia- 
tion. The FEL collective instability in the electron beam produces an exponential growth (along 
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the undulator) of the modulation of the electron density on the scale of undulator radiation 
wavelength. The fluctuations of current density in the electron beam are uncorrelated not only 
in time but in space, too. Thus, a large number of transverse radiation modes are excited when 
the electron beam enters the undulator. These radiation modes have different gain. As undulator 
length progresses, the high gain modes start to predominate more and more. For enough long 
undulator, the emission will emerge in a high degree of transverse coherence. An intensity gain 
in excess of 10^ — 10^ is obtained in the saturation regime. At this level, the shot noise of the 
electron beam is amplified up to complete micro-bunching, and all electrons radiate almost in 
phase producing powerful, coherent radiation. 

Understanding of coherence properties of the radiation from SASE FEL is of great practical im- 
portance. Properties of the longitudinal coherence have been studied in [12,13,14,15,16,17,18]. 
It has been found that the coherence time increases first, reaches maximum value in the end of 
the linear high gain regime and then drops when amplification process enters nonlinear stage 
[18]. The first analysis of the problem of transverse coherence has been performed in [19]. The 
problem of start-up from the shot noise has been studied analytically and numerically for the 
linear stage of amplification. It has been found that the process of formation of transverse coher- 
ence is more complicated than that given by naive physical picture of transverse mode selection. 
Namely, even after finishing the transverse mode selection process the degree of transverse co- 
herence of the radiation from SASE FEL visibly differs from the unity. This is consequence 
of the interdependence of the longitudinal and transverse coherence. The SASE FEL has poor 
longitudinal coherence which develops slowly with the undulator length thus preventing a full 
transverse coherence. First studies of the evolution of transverse coherence in the nonlinear 
regime of SASE FEL operation have been performed in [20]. It has been found that similarly 
to the coherence time, the degree of transverse coherence reaches maximum value in the end of 
the linear regime. Further increase of the undulator length leads to its decrease. Despite output 
power of the SASE FEL grows continuously in the nonlinear regime, maximum brilliance of 
the radiation is achieved in the very beginning of the nonlinear regime. Due to a lack of com- 
puting power available at that time we limited our study with a specific numerical example just 
illustrating the general features of coherence properties of the radiation produced by the SASE 
FEL operating in the nonlinear regime. 

In this paper we present general analysis of the coherence properties (longitudinal and trans- 
verse) of the radiation from SASE FEL. The results have been obtained with time-dependent, 
three-dimensional FEL simulation code FAST [21] performing simulation of the FEL process 
with actual number of electrons in the beam. Using similarity techniques we present universal 
dependencies for the main characteristics of the SASE FEL covering all practical range of X-ray 
FELs. 
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2 Basic relations 



Design of the focusing system of XFEL assumes nearly uniform focusing of the electron beam 
in the undulator, so we consider axisymmetric model of the electron beam. It is assumed that 
transverse distribution function of the electron beam is Gaussian, so rms transverse size of 
matched beam is a = y/ej3 ,where e = en/7 is rms beam emittance and P is focusing beta- 
function. An important feature of the parameter space of XFEL is that the space charge field 
does not influence significantly on the FEL process and calculation of the EEL process can be 
performed by taking into account diffraction effects, the energy spread in the electron beam, 
and effects of betatron motion only. In the framework of the three-dimensional theory operation 
of the EEL amplifier is described by the following parameters: the diffraction parameter B, the 
energy spread parameter A^, and the betatron motion parameter A;^ [22,23]: 



kl^{a^/Ef/p\ (1) 



where V = [Iuj'^6'^A'j;j/{lAc'^'-f'^'-f)] is the gain parameter and p = cj'^T/uj is the efficiency 
parameter ^ . When describing shot noise in the electron beam, one more parameter appears, 
the number of electrons on the coherence length, = I / {eujp). The following notations are 
used here: / is the beam current, uj = 2nc/X is the frequency of the electromagnetic wave, 
= -ft'rms IS the mis uudulator parameter, 7 is relativistic factor, 7^^ = 7^^ + 9'^, 

= 27r/A„ is the undulator wavenumber, /a = 17 kA is the Alfven current, Ajj = 1 for 
helical undulator and Ajj = Jo{K'^j^^^/2{l+K'^^g)) — Ji{K'^^^g/2{l+K'^^^^^)) for planar undulator. 
Here Jo and Ji are the Bessel functions of the first kind. The energy spread is assumed to be 
Gaussian with rms deviation ue. 

The amplification process in the FEL amplifier passes two stages, linear and nonlinear. The 
linear stage lasts over significant fraction of the undulator length (about 80%), and the main 
target for XFEL optimization is the field gain length. In the linear high-gain limit the radiation 
emitted by the electron beam in the undulator can be represented as a set of modes: 

+ iEy^ / do; ex.p[iu;(z/c - t)] x ^nm(<^, z)^nm('r, u) exp[A„„(a;)2; + incj)] . (2) 

n,m 

^ Note that it differs from 1-D definition by the factor B^/^ [22]. 
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When amplification takes place, the mode configuration in the transverse plane remains un- 
changed while the amplitude grows exponentially with the undulator length. Each mode is char- 
acterized by the eigenvalue A„^(a;) and the field distribution eigenfunction ^nm{r, co) in terms 
of transverse coordinates. At sufficient undulator length fundamental TEMqo mode begins to 
give main contribution to the total radiation power. Thus, relevant value of interest for XFEL 
optimization is the field gain length of the fundamental mode, Lg — 1/ Re(Aoo), which gives 
good estimate for expected length of the undulator needed to reach saturation, Lgat — 10 x Lg. 
Optimization of the field gain length is performed by means of numerical solution of the corre- 
sponding eigenvalue equations taking into account all the effects (diffraction, energy spread and 
emittance) [23,24]. Computational possibilities of modem computers allows to trace complete 
parameter space of XFEL (which in fact is 11 -dimensional). From practical point of view it is 
important to find an absolute minimum of the gain length corresponding to optimum focusing 
beta function. For this practically important case the solution of the eigenvalue equation for 
the field gain length of the fundamental mode and optimum beta function are rather accurately 
approximated by [25]: 



/ T'. \ -.3/2 A 1/2 



where a^ — a^j m^c? . Accuracy of this fit is better than 5 % in the range of parameter e = 27re/ A 
from 1 to 5. 

Equation (3) demonstrates clear interdependence of physical parameters defining operation of 
the XFEL. Let us consider the case of negligibly small energy spread. Under this condition 
diffraction parameter B and parameter of betatron oscillations, A;^ are functions of the only 
parameter e: 

S~12.5xe^/% = l/(/3r) ~ 0.158 X e^/^ (4) 



EEL equations written down in the dimensionless form involve an additional parameter A^'c 
defining the initial conditions for the start-up from the shot noise. Note that the dependence of 
output characteristics of the SASE EEL operating in saturation is slow, in fact logarithmic in 
terms of N^- Thus, we can conclude that with logarithmic accuracy in terms of N^, character- 
istics of the SASE EEL written down in a normalized form are functions of the only parameter 
e. 
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3 An approach for numerical simulations 



Rigorous studies of the nonlinear stage of amplification is possible only with numerical simula- 
tion code. Typically FEL codes use an artificial ensemble of macroparticles for simulating the 
FEL process when one macroparticle represents large number of real electrons. Thus, a natural 
question arises if macroparticle phase space distributions are identical to those of actual electron 
beam at all stages of amplification. Let us trace typical procedure for preparation of an artifi- 
cial ensemble [26,27]. The first step of particle loading consists in a quasi-uniform distribution 
of the macroparticles in the phase space. At this stage an ensemble of particles with random 
distribution is generated which occupies a fraction of the phase space. Then this ensemble is 
copied on the other parts of the phase space to provide pseudo-uniform loading of the phase 
space. Pseudo-uniformity means that initial microbunching at the fundamental harmonic (or for 
several harmonics) is equal to zero. Also, phase positions of the mirrored particles are corre- 
lated such that microbunching does not appear due to betatron oscillations, or due to the energy 
spread. Finally, artificial displacements of the macroparticles are applied to provide desired (in 
our case gaussian) statistics of microbunching at the undulator entrance. We note that it is not 
evident that such an artificial ensemble reflects actual physical situation for a short wavelength 
SASE FEL. Let us consider an example of the SASE FEL operating at the radiation wavelength 
of 0.1 nm. With the peak current of 5 kA we find that the number of electrons per wavelength 
is about 10^. On the other hand, it is well known that properties of an artificial ensemble (even 
at the first step of pseudo-uniform loading) converge very slowly to the model of continuous 
media. In fact, even with the number of macroparticles per radiation wavelength 6.4 x 10"^ the 
FEL gain still visibly deviates from the target value. Introducing of an artificial noise makes sit- 
uation with the quality of an ensemble preparation even more problematic. The only way to test 
the quality of an artificial ensemble is to perform numerical simulations with actual number of 
electrons in the beam. We constructed such a version of three-dimensional, time-dependent FEL 
simulation code FAST [21]. Comparison of the results with direct simulations of the electron 
beam and with an artificial distributions has shown that artificial ensembles are not adequate to 
the problem. Artificial effects are pronouncing especially when calculating such fine features 
as transverse correlation functions. Thus, all the simulations presented in this paper have been 
performed with code FAST using actual number of electrons in the beam. 

4 General overview of the properties of the radiation form SASE FEL 

The result of each simulation run contains an array of complex amplitudes E for electromag- 
netic fields on a three-dimensional mesh. At the next stage of the numerical experiment the data 
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arrays are handled with postprocessor codes to calculate different characteristics of the radia- 
tion. However, as the first step it is worthwhile to obtain qualitative analysis of the object under 
study. The plots in an upper row of Fig. 1 show evolution of the power density distribution, 
/ = in a slice of the radiation pulse along the undulator. We see that due to the start-up 
of amplification process from the shot noise many transverse radiation modes are excited when 
electron beam enters the undulator. Mode selection process (2) serves as a filter for selection of 
the fundamental radiation mode having maximum gain. 

Integration of the power density over transverse cross section of the photon beam gives us 
instantaneous radiation power, P oc J I dr±- Evolution of temporal structure of the radia- 
tion power along the undulator is traced in Fig. 1 in terms of normalized radiation power 
77 = P/{pPb) where Pb — ^mc^I/e is the electron beam power. Averaging of the radiation 
power along the pulse gives us averaged radiation power. Evolution of normalized averaged 
power (77) along normalized undulator length z — Tzis shown in Fig. 2. Note that the radiation 
produced by SASE FEL operating in the linear regime holds properties of completely chaotic 
polarized light [18] - a statistical object well described in the framework of statistical optics 
[28]. This is simple consequence of the fact that the shot noise in the electron beam is a Gaus- 
sian random process. The FEL amplifier, operating in the linear regime, can be considered as a 
linear filter which does not change the statistics of the signal. As a result, we can define general 
statistical properties of the output radiation without any calculations. For instance, in the case 
of the SASE FEL the real and imaginary parts of the slowly varying complex amplitudes of the 
electric field of the electromagnetic wave, E , have a Gaussian distribution, the instantaneous 
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Fig. 1. Evolution of the power density distribution in a slice (top) and normalized power in a radiation 
pulse (bottom) versus s = pooq{z/vz — t) for the reduced lengths z = 18.4, 31.6, and 44.8 (left, middle, 
and right plots, respectively). Here e = 2. Crosses show geometrical center of the radiation beam. 
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Fig. 2. Averaged normalized efficiency, ?}, and normalized rms deviation of instantaneous radiation 
power, (Tp along the normalized undulator length, z = Fz. Here e = 2. 




Fig. 3. Right plot: normaUzed coherence time fc, along the normalized undulator length, z = Tz. 
Left plot: degree of transverse coherence, ( (soUd Une), and normalized degeneracy parameter of the 
radiation, S (dashed line), along the normaUzed undulator length, Tz. Dotted line shows Up (see Fig. 2). 
Here e = 2. 

power density, / = fluctuates in accordance with the negative exponential distribution 
(see Fig. 4): 



Due to the start-up of amplification process from the shot noise many transverse radiation modes 
are excited when electron beam enters the undulator. For gaussian random process any integral 
of the power density, for example, radiation power P, fluctuates in accordance with the gamma 
distribution: 

/ P \ ^-^ 1 ( P \ 
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Fig. 4. Probability density distributions of the instantaneous power density / = |£'p (top), and of the 
instantaneous radiation power P (bottom) from SASE FEL for different stages of amplification: linear 
regime (z = 32), saturation regime (z = 63), and deep nonlinear regime (z = 105). Solid lines on 
the power density histograms (top) represent negative exponential distribution (5). Solid lines on power 
histograms (bottom) represent gamma distribution (6) with M = l/ap (see Fig. 2). Here e = 2. 

where r(M) is the gamma function with argument M, and 

M=\, (7) 



and dp = ((P — (P))^)/(P)^ is the relative dispersion of the radiation power. Note that for 
completely chaotic polarized light parameter M has clear physical interpretation, it is the num- 
ber of modes [22]. Thus, it becomes clear that the relative dispersion of the radiation power 
directly relates to the coherence properties of the SASE FEL operating in the linear regime. The 
degree of transverse coherence in this case can be naturally defined as: 

( = -L = ^2 _ (g) 



Indeed, in the linear regime we deal with a Gaussian random process, and the power density 
fluctuates in accordance with the negative exponential distribution and its relative width is equal 
to 1 . If there is full transverse coherence then the same refers to the power. If the radiation is 
partially coherent, then we have a more general law for power fluctuations, namely the gamma 
distribution (6). In the linear regime fluctuations of the radiation power, crp, grow steadily with 
the undulator length (see Fig. 2) because of mode selection process (2). 

Another physical parameter of the problem relating to the transverse coherence is diffraction pa- 
rameter B. Mechanism of formation of transverse coherence is rather transparent. If diffraction 
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expansion of the radiation on a scale of the field gain length is comparable with the transverse 
size of the electron beam, we can expect a high degree of transverse coherence. For this range 
of parameters the value of the diffraction parameter is small. If diffraction expansion of the 
radiation is small (which happens at large values of the diffraction parameter) then we can ex- 
pect significant degradation of the degree of transverse coherence. This happens simply because 
different parts of the beam produce radiation nearly independently. In terms of the radiation ex- 
pansion in the eigenmodes (2) this range of parameters corresponds to degeneration of modes 
(see Appendix A). Diffraction parameter for optimized XFEL exhibits strong dependence on the 
parameter e (see eq. (4) and Fig. 5), and we can expect that the degree of transverse coherence 
should drop rapidly with the increase of the parameter e. 

An important physical quantity describing random fields is the coherence time. Figure 1 gives 
qualitative picture of formation of longitudinal coherence in SASE FEL. At the beginning of 
the undulator the radiation is simply incoherent undulator radiation. When amplification process 
starts to dominate over spontaneous emission, we obtain formation of spikes (wavepackets). The 
width of the spikes defines the coherence time. In the high gain linear regime the width of spikes 
grows with the undulator length, and the coherence time also grows proportionally to the square 
root of the undulator length. It achieves maximum value in the end of linear regime and then 
decreases rapidly in the nonlinear regime (see Fig. 3). The maximum value of the coherence 
time depends on the saturation length and, therefore, on the value of the parameter Nc [18]. 

We see that physical background defining general features of the radiation from SASE FEL 
operating in the linear regime is rather transparent. Despite the behavior of the SASE FEL in 
the nonlinear regime is rather complicated (see Figs. 2 and 3), we will show below that the main 
characteristics of the SASE FEL operating at the saturation point have rather simple physical 
scaling. 
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Fig. 5. Diffraction parameter B versus parameter e for optimized XFEL. 
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5 General definitions 



5.1 Degree of transverse coherence 



The transverse coherence properties of the radiation are described in terms of the transverse 
correlation functions. The first-order transverse correlation function is defined as 



{E{f^,z,t)E*{ff^,z,t)) 
{\E{f^,z,m{\E{ff^,z,t)\^f^' 



where E is the slowly varying amplitude of the amplified wave: 

E = E{r^, z, t)e'^o{z/c-t) ^ q q _ 



In the following we consider the model of a stationary random process, thus assuming that 71 
does not depend on time. We define the degree of transverse coherence as: 



When SASE PEL operates in the linear regime such a definition for the degree of transverse 
coherence is mathematically equivalent to (8) expressed in terms of the relative dispersion of the 
instantaneous radiation power, (Xp. Analysis of the asymptotic of the deep nonlinear regime (see 
Fig. 3) shows that surprisingly the degree of transverse coherence defined by (10) again tends 
to be an agreement with (8). This feature has deep physical background. When amplification 
process just enters nonlinear stage, the statistics of the radiation is not Gaussian anymore. In 
particular, the probability distribution function of the radiation power density, is not negative 
exponential distribution (see Fig. 4). Thus, definition of the degree of transverse coherence (8) 
has no physical sense near the saturation point. However, simulations show that in the deep 
nonlinear regime the probability distribution of the radiation power density again tends to the 
negative exponential distribution. This gives us a hint that the properties of the radiation from 
SASE EEL operating in the deep nonlinear regime tend again to be the properties of completely 
chaotic polarized light. Similar asymptotical behavior has been observed in the framework of 
the one-dimensional model as well [18]. 
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5.2 Coherence time 



Longitudinal coherence is described in terms of time correlation functions. The first order time 
correlation function, gi{t, t'), is calculated in accordance with the definition: 



gi{r.t-t') 



{E{r,t)E*{f,t')) 



{\E{f,t) \'){\mt') 



1/2 



(11) 



For a stationary random process time correlation functions are functions of the only argument, 
T = t — t'. The coherence time is defined as [22]: 



oo 

/ 



l^iWI dr. 



(12) 



5.3 Brilliance and degeneracy parameter 



Main figure of merit of the SASE FEL performance is brilliance, i.e. density of photons in the 
phase space. In fact, the brilliance is proportional to the degeneracy parameter 5, i.e. the number 
of photons per mode (coherent state). Note that when 5^1, the classical statistics is applicable, 
while quantum description of the field is necessary as soon as 5 is comparable to (or less than) 
one. Using the definitions of the degree of transverse coherence (10) and coherence time (12), 
one can define degeneracy parameter: 

8 = NphTcC , (13) 



where Nph — N^^^/T is the photon flux, N^'f^ is the total number of photons in a long flat- 
top pulse of duration T (as everywhere in this paper we consider ensemble average values). 
The definition (13) is perfect for a Gaussian random process (that we have in linear regime). 
Indeed, degree of transverse coherence is an inverse number of transverse modes (8), while 
Tc/T is an inverse number of longitudinal modes within the pulse [22,28]. Thus, degeneracy 
parameter is equal to the number of photons per pulse divided by total number of modes per 
pulse (that is equal to the squared inverse rms fluctuations of pulse energy [22]). It can be 
directly measured in an experiment for SASE FEL operating in the linear regime. We use (13) 
as a natural generalization to characterize SASE FEL properties at saturation. 

Let us introduce a notion of normalized degeneracy parameter 

6 = fiU . (14) 
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Here normalized FEL efficiency is defined as r) — P/ (pPh) where P is radiation power, and 
Pb = ^mc^I/e is electron beam power. Normalized coherence time is defined as fc = pcoTc. 
Parameter S and the degeneracy parameter 5 are simply related as: 

' 5 , (15) 



or, in practical units 

5 ~ 2.7 X 10^ X A^[A ] X /[kA] X £;[GeV] x 5 , (16) 

where E — ^mc^ is the electron energy. Note that the degeneracy parameter is very large even 
for a S ASE FEL operating at the wavelength of 0. 1 nm. With multi-kA electron beams and other 
relevant parameters, mentioned in Table 1, degeneracy parameter would be of the order of 10^ 
- 10^. Thus, a classical treatment of SASE FEL is justified. 

Let us now turn to the calculation of peak brilliance. It is defined as a transversely coherent 
spectral flux: 

^^^^ (17) 

The spectrum of SASE FEL radiation in a high-gain linear regime has a Gaussian shape, it is 
also close to the Gaussian at saturation point [22]. In this case 



where a,^ is the rms bandwidth. For a Gaussian line, with the definition of coherence time (12), 
one gets [22]: 

'c — 

Thus, it follows from (14) and (17) that 

B, = i^^. (18) 

The peak brilliance can then be calculated as follows 

S,[phot./(sec. mrad" mm' 0.1% bandw.)] ~ 4.5 x 10^^ x ^ ■^[^'''^^ ^ ^ (19) 

A[A J 
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For future SASE FELs, operating at 1 A , an expected value of the peak brilliance is 10^^-10^^. 

An important feature of our analysis is application of similarity techniques to the analysis of 
simulation results which allows to derive universal parametric dependencies of the output char- 
acteristics of the radiation. As we mentioned above, within accepted approximations (optimized 
SASE PEL and negligibly small energy spread in the electron beam), output characteristics of 
SASE PEL are universal functions of two parameters, normalized undulator length z — Tz and 
parameter e. If one traces evolution of the brilliance (or, degeneracy parameter) of the radiation 
along the undulator length there is always the point, which we define as the saturation point, 
where the brilliance reaches maximum value (see Fig. 3). In the following we present character- 
istics of the radiation at the saturation point which are universal functions of the only parameter 
e. 

6 Properties of the radiation from optimized XFEL operating in the saturation regime 

Simulations of the FEL process have been performed for the case of a long bunch with uniform 
axial profile of the beam current. Such a model provides rather accurate predictions for the 
coherence properties of the XFEL, since typical radiation pulse from the XFEL is much longer 
than the coherence time. Calculations has been performed with FEL simulation code FAST 
using actual number of electrons in the beam. The value of parameter A'c = 8 x 10^ corresponds 
to the parameter range of XFEL operating at the radiation wavelength about 0.1 nm. 




Fig. 6. Top: normalized power in a radiation pulse, fj versus s = pu{z/vz — t). Bottom: typical power 
density distribution in a slice of the radiation. SASE FEL operates in the saturation. Values of e = 1, 2 
and 4.5 correspond to the left, middle, and right plot, respectively. Crosses show geometrical center of 
the beam. 
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Fig. 7. Saturation length Zgat = Tzgat versus parameter e. 
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Fig. 8. Degree of transverse coherence, Qax^ and normaUzed coherence time, in the saturation versus 
parameter e. 
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Fig. 9. Averaged efficiency, (7}sat)> and normahzed degeneracy parameter, ^gat, in the saturation versus 
parameter e. 

Figure 6 gives visual picture of the slice properties of the radiation at the saturation for different 
values of the parameter e. Saturation point is defined as the point where the radiation achieves 
maximum brilliance (or, maximum degeneracy parameter (14)). A series of simulation runs 
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Fig. 10. Partial contributions of the modes with azimuthal index m = . . . 4 into the total power versus 
parameter e. SASE FEL operates in the saturation. 

has been performed in the range of the parameter e = 0.25. ..4.5. Application of similarity 
techniques described above allowed us to extract universal parametric dependencies of the main 
characteristics of the optimized XFEL operating in the saturation regime (see Figs. 7-9). 

Figure 7 shows the dependence of the saturation length ^sat = Tzsat on parameter e. Analysis of 
the curve shows that the saturation length scales as 5sat oc e^/^. Such dependence directly fol- 
lows from the optimization procedure of the gain length given by (3). The normalized coherence 
time in the saturation regime, f^*^* is also proportional to e^/^ (see Fig. 8). 

The dependence of the degree of transverse coherence in the saturation regime on the parameter 
e exhibits rather complicated behavior (see Fig. 8). It reaches maximum value in the range of e 
values about of unity, and drops at small and large values of e. Actually, the degree of transverse 
coherence is formed due to two effects. The first effect takes place due to interdependence of the 
poor longitudinal coherence and transverse coherence [19]. Due to the start-up from shot noise 
every radiation mode entering (2) is excited within finite spectral bandwidth. This means that 
in the high gain linear regime the radiation of the SASE FEL is formed by many fundamental 
TEMoo modes with different frequencies. The transverse distribution of the radiation field of the 
mode is also different for different frequencies. Smaller value of the diffraction parameter (i.e. 
smaller value of e) corresponds to larger deviation of the radiation mode from the plane wave. 
This explains a decrease of the transverse coherence at small values of e. When the parameter 
e increases, the diffraction parameter increases as well thus leading to the degeneration of the 
radiation modes. Amplification process in the SASE FEL passes limited number of the field 
gain lengths, and starting from some value of e the linear stage of amplification becomes too 
short to provide mode selection process (2). When amplification process enters nonlinear stage, 
the mode content of the radiation becomes even more rich due to independent growth of the 
radiation modes in the nonlinear medium (see Fig. 10). Thus, at large values of e the degree 
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Table 1 



Parameter space of XFEL projects 





European XFEL 


LCLS 


scss 




SASEl 


SASE2 






Radiation wavelength, nm 


0.1 


0.15 


0.15 


0.1 


Beam energy, GeV 


17.5 


17.5 


14.35 


6.135 


rms normalized emittance e„, mm-mrad 


1.4 


1.4 


1.2 


0.85 


Parameter e 


2.6 


1.7 


1.8 


4.5 


Degree of transverse coherence ^ 


0.65 


0.85 


0.83 


0.24 



of transverse coherence is limited by poor mode selection. Analytical estimations, presented in 
Appendix A, show that in the limit of large emittance, e » 1, the degree of transverse coherence 
scales as 

We present in Fig. 9 the plots for normalized efficiency and degeneracy parameter for optimized 
XFEL. Normalized efficiency in saturation has simple scaling, it falls inversely proportional to 
the parameter e. Taking into account that the value of the coherence time f^^* scales proportional 
to e^/^, we find that the normalized degeneracy parameter of the radiation is nearly proportional 
to the degree of transverse coherence, Sgat oc C/^^^^- 

Finally, in Table 1 we present comparison of existing XFEL projects, the European XFEL, 
LCLS and SCSS in terms of degree of transverse coherence [4,5,6]. We see that the European 
XFEL and LCLS are in the same range of parameter space. These projects assume conserva- 
tive value of the emittance, and relatively high degree of transverse coherence is achieved by 
increasing the energy of the driving accelerator. Project SCSS assumes much smaller energy of 
the driving accelerator. Thus, despite much smaller value of the normalized emittance it falls in 
the range of parameters for the output radiation with poor transverse coherence. 
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A Solution of the eigenvalue equation and estimates of transverse coherence in the limit 
of wide electron beam 



A.1 Basic equation 



Let us have at the undulator entrance a continuous electron beam with the current /q, with the 
Gaussian distribution in energy 



(A.1) 



and in a transverse phase plane 



+ {x'f/k} 

2^2 



(A.2) 



the same in y phase plane. Here kjj — 1//3 is the wavenumber of betatron oscillations and 

The results of this paper can be used in the case of superposition of the natural undulator focus- 
ing and an alternating-gradient external focusing if the following condition is satisfied [23]: 

<^ mm 1, ' 



2Tre ' 



where L/ is a period of the external focusing structure, /3 is an average beta-function, e is rms 
emittance of an electron beam, and A is a radiation wavelength. This condition is met in many 
practical situations. 

Using cylindrical coordinates, in the high-gain limit we seek the solution for a slowly varying 
complex amplitude of the electric field of the electromagnetic wave in the form [22]: 

(sm(rnp) \ 
, (A.3) 
cos(n</?) / 



where n is an integer, n > 0. For each n there are many radial eigenmodes that differ by 
eigenvalue A and eigenfunction $„„j(r). The integro-differential equation for radiation field 
eigenmodes [24,31,32] taking into account the space charge effect [23] can be written in the 
following normalized form: 
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-4 J df'f' {^r.m{r') 



til 

2 



+ 



1 d 



n 



df'2 ' r'df'"f^ + ^^^'^ 



X 



/dC. 

SI 



c 



sin (kpC) 



exp 



-(A + iC)C 



X exp 



x/„ 



(l-i^fe^C/2)(r^ + r^^) 
sin^(^/3C) 

2(1 -iS^^C/2)rf'cos(^^C) 
sin^(A;/3C) 



(A.4) 



where /„ is the modified Bessel function of the first kind. The following notations are used 
here: f = r/{a\^), B = 2a'^Vijj/c is the diffraction parameter, = kp/V is the betatron 
motion parameter, Ap = 2c^(Ajj^^sO'^^)"^ is the space charge parameter, A| = ((A£^)^) / {p^S"^) 
is the energy spread parameter, C = [k^ — a;/(2c7^)] /F is the detuning parameter, F = 

1 /2 

AjjIoLu'^d'^ (/ac^7^7)~ is the gain parameter, p = cy'^T/uj is the efficiency parameter, uj is 
the frequency of the electromagnetic wave, = ^rms/ 7, ^rms is the rms undulator parameter, 7 
is relativistic factor, 7"^ = 7"^ + , /c„ is the undulator wavenumber, /a = 17 kA is the Alfven 
current, Ajj = 1 for helical undulator and Ajj = Jo{K^^j2{l+K^^J)-Ji{K^^j2{l+K^^J) 
for planar undulator. Here Jq and Ji are the Bessel functions of the first kind. The space charge 
effect is included into (A.4) under the condition cr^ » 0^7^/0;^. 

A.2 Exact solution 

As suggested in [24] we apply to (A.4) the Hankel transformation defined by the following 
transform pair: 

00 00 

^nm(p) = j dffJn{pf)^nm{r) , ^nm(^) = j dppJn{pr)^nm{p) ■ 



Then we obtain the integral equation for the Hankel transform ^nm{p) [23]: 



*nm(p) 



2 i SA - p2 



dp'p'^nmip') 



X 





„ (i-iBk}c/2y 



^ ^ Al{2iBA-p'') 



c 



exp 



A^C^ 



(A + iC-X 
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X exp 



p2 _^ 



4(l-i5HC/2) 



pp' cos(A;gC) 
2(l-ii?HC/2) 



(A.5) 



When the space charge field is negligible, Ap — > 0, this equation is reduced to that derived in 
[24]. 

To solve (A.5) we discretize it: 

p'.^Aij-^), j^l,2,...,K, 

where A and K should be chosen such that the required accuracy is provided. Then we obtain 
a matrix equation 



or, [Mn — I]^nm — 0, wherc / is a unit matrix. Matrix M„ depends on an eigenvalue A as 
well as on the problem parameters: B, k^, A|, Ap, and C. The eigenvalues of all radial modes 
for a given azimuthal index n can be found by solving the equation |M„ — /| = 0. Then the 
calculation of the eigenmodes is straightforward. 

This algorithm allows one to find with any desirable accuracy the eigenvalues and eigenfunc- 
tions of a high-gain FEL including all the important effects: diffraction, betatron motion, energy 
spread, space charge, and frequency detuning. Therefore, it can be considered as a universal tool 
for calculation and optimization of high-gain FELs of wavelength range from infrared to X-ray. 



A.3 Parallel cold beam, large diffraction parameter 



A parallel beam limit is the case when there is no betatron motion, i.e. kp — > 0. Let us also 
assume here for the sake of compactness that the effects of energy spread and space charge can 
be neglected (A^ — > 0, Ap — > 0). Equation (A.4) can then be reduced to 

To find explicit solutions of the eigenvalue equation in the limit of large diffraction parameter 
(more specifically, we require B^l^ ^ 1), we use the variational method [24,29]. We construct 
a variational functional from (A.6): 



^nmif) = -^-^^^nm{r) . (A.6) 
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j d 



nm\1") 



dr^ fdr 

;2 

(A + iC) 



-/drf^l^^$L(r), (A.7) 



and seek for a solution in the form 

^nm{r) = r'exp{-af^)Ll{2af^) , (A.8) 

where are associated Laguerre polynomials 

ml fc! Vm — k/ 

and 

'n\ n! 



/ty k\{n-k)\ 



is a binomial coefficient. 

Equation (A.7) and the variational condition, 6 A/ 5a = 0, lead to the following two equations 
for two unknown quantities, A and a: 

iBA 1 / l + /; + 2m\ 

l + n + 2m-'- 1 =0, (A.10) 

a a(A + iC)2 V 2a J 

i5A+,_l^ri-i±^Uo. (A.11) 



(A + iC)^ 

Solving Eqs. (A. 10), (A. 11) in zero order, we get 1-D asymptote for the eigenvalue equation 
[22]. The eigenvalue for a growing solution reaches maximum at (7 = 0: 

Then we can find first order correction (in B^^^^) to the growth rate Re A, and a mode parameter 



25V3 I 3^3^1/3 i ' ^ ^ 



20 



2^2 



(A. 13) 



Equations (A. 8), (A. 12), and (A. 13) are the solutions for field distributions and growth rates of 
eigenmodes of a high-gain PEL with a cold parallel beam in the limit B^/^ » 1. Note that in 
[22,30] the exact solution of the eigenvalue problem for a parabolic beam density distribution 
was obtained. The eigenfunctions were expressed in terms of the confluent hypergeometric 
function. In this case, in the limit of a large diffraction parameter the growth rates of eigenmodes 
are reduced to (A. 12). The confluent hypergeometric function takes the form of the associated 
Laguerre polynomials, so that field distribution is reduced to (A. 8) with the parameter a given 
in (A. 13). This is not a surprise because in this limit the width of the field distribution is much 
smaller than the width of the electron beam (Re a ~ B^/^ » 1), and the electron density 
function behavior is important only near the axis. This behavior is the same (quadratic) for 
parabolic and Gaussian distributions. We can also conclude that this asymptotical solution is 
valid for any axisymmetric density distribution with the maximum density on axis. We also 
see that the variational solution is a good asymptotical method. The attempts to generalize it to 
the entire (or, at least wider) parameter range [24,29,33] lead to the loss of accuracy control, 
although can give a practically useful fit of the exact solution within some range. 



A.4 Large emittance 



Let us still assume that the space charge effect can be neglected (Ap 0). Applying now the 
variational method to the Eq. (A.4) with the trial functions (A.8), we obtain for large diffraction 
parameter: 



\BK If, 

1 + n + 2m / dC 

a a J 1 — i 



c 



X 



1 - 



(1 + 71 + 2to) 



SHC/2 



-\Bk%/2 



exp 



+ 



yBk}C/2 



0, 



(A. 14) 



iBK + 



oo 



_c 



exp 



-(A + iC'K 



X 



'l + n + 2m){l-iBklC/2) 



. 



(A. 15) 



The system of equations (A. 14) and (A. 15) can be solved numerically. In the following we ne- 
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gleet the effect of the energy spread, — > 0. We also assume that beta-function is optimized 
for the highest FEL gain as it happens in practice. Since diffraction parameter depends on beta- 
function, it is more convenient go over to other normalized parameters. Indeed, diffraction pa- 
rameter can be rewritten as B — 2e//c/j, where e = 27re/A. Then we can go from parameters 
{B, kp) to (e, fc/j). After optimizing parameters kp and C, we will find growth rates and eigen- 
functions of all eigenmodes depending on the only parameter, e. Equations (A. 14) and (A. 15) 
can now be written as (e ^ 1): 



1 + n + 2m 



2ieA 

akf- 



X 



'0 
n-\- 2m) 



a J 1 - i 



ekpC 



exp 



-(A + i(5)C) 



ak^e 



iekpC 



+ 



1 - i ekpC 



= 0, 



(A.16) 



1 - i ekpC, 



exp 



(l + n + 2m)(l-iefeX) ^ 



In zero order we find [24]: 



= 0. (A. 17) 



An 



(At ^ 



exp 



-(Ao + iC^)C 



(A.18) 



Solving this equation numerically, we find that maximal growth rate 

T. A 0-37 
ReAo ~ —:— 



is achieved at the optimal values of kp ~ 0.5/e^ and C ^ 0.4/e. Note that for optimal beta- 
function the diffraction parameter can be expressed as S ~ 4e^. 



Solving Eqs. (A.16) and (A. 17) in the first order in e ^, we obtain: 



0.37 / 0.83(1 + 71 + 2m) ^ 



a~ (0.44- 0.51 i)e. 



(A. 19) 
(A.20) 



Equations (A. 8), (A. 19), and (A.20) are the solutions for the field distributions and growth rates 
of eigenmodes of a high-gain FEL with optimized beta-function in the limit of e 1. 



22 



A. 5 Estimates of transverse coherence 



We can now make a simple estimate of the number of transverse modes, M, in the high-gain 
linear regime of a SASE PEL operation. The degree of transverse coherence would then be the 
inverse number of modes (see (8)): 

C = = 4 , (A.21) 

where ap = ((P — {P)Y) / {P^ is the relative dispersion of the radiation power. The field of 
the electromagnetic wave can be represented as a set of modes, see (2). In the limit, considered 
here {B^/^ ^ 1 or e 1), the modes are orthogonal, and the total power can be written as ^ 

OC' oo 

Ptot{z) ^Y.Pnm{z) = 2J2Y1 ^nm exp(2 Re A^^f) - Aqo exp(2 Re Aooi) • (A.22) 

Summation over azimuthal index n is done twice here since for n 7^ there are two orthogonal 
modes that degenerate [22]. One can also show that (here we assume it without a proof) the 
factors Anm are the same for all modes in the considered asymptote ^ . Since the power of each 
mode fluctuates in accordance with the negative exponential distribution (5), the dispersion is 
equal to the squared average power for each mode. The total dispersion is simply the sum of 
dispersions because the modes are independently excited. Thus, the inverse relative dispersion 
(or, number of modes) can be calculated as (E PumT/ Yl Pnm' explicitly: 

/ 00 00 \ 2 

2 E E exp[-2Ngb{n + 2m)] - 1 

M ^ -=°-=" ^ , (A.23) 

2 E E exp[-4:Ngb{n + 2m)] - 1 

n=0 m=0 



where Ng — Re Aoi is a number of field gain lengths along the undulator, b — \f2l (3V^S^/^) 
for a cold parallel beam, and b — 0.83/e for a beam with large emittance and optimized beta- 
function. Equation (A.23) is valid when 6 <^ 1 and Ng 1. 

In a particular case when 1 ^ A^^ ^ 6^ the summation in (A.23) can be substituted by the 
integration. Then for a cold parallel beam we get: 

27 //?^/'^\^ 

M ~ — for K iV„ < B^/^ . (A.24) 

2 \Ng ) 

^ Ensemble average is meant here. 

^ More strictly, orthogonality of the radial modes with the same azimuthal index, as well as equality of 
the factors Anm, hold with an accuracy e^^ ^ 1. Taking these corrections into account would result in 
the correction of the order of (Re Aoi)~^ to the number of modes. 
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For a beam with a large emittance and optimized beta-function the number of modes is 



We note that applicability region of these estimations is the high-gain linear regime. Numerical 
simulations presented in this paper show that the maximum degree of transverse coherence is 
achieved already in the nonlinear mode of operation. Linear analysis, presented here, does not 
allow to describe this maximum degree of transverse coherence. However, it can be roughly 
estimated if one substitutes Ng by the number of field gain lengths at the end of the linear 
regime. As an estimate, one can take about 70% of the number of field gain lengths required to 
reach saturation . In any case the asymptotical behavior of the degree of transverse coherence 
is 



in the case of a beam with large emittance and optimized beta-function. 



^ Note that saturation occurs earlier for a larger number of modes. This would give a weak (logarithmic) 
correction to the value of the transverse coherence. 




(A.25) 
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